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Abstract 

Jet finding algorithms, as they are used in e~^e~ and hadron collisions, are 

reviewed and compared. It is suggested that a successive combination style 

algorithm, similar to that used in e~^e~ physics, might be useful also in hadron 

collisions, where cone style algorithms have been used previously. 
12.38.Bx, 13.87.-a 
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I. INTRODUCTION 



The measurement of jet cross sections has provided useful tests of Quantum Chromody- 
namics both at hadron coUiders and at electron-positron colliders. The observed jets provide 
a view of the underlying hard quark and gluon interactions that occur at very small distance 
scales. However, this view is inevitably clouded by the subsequent long distance shower- 
ing and eventual hadronization of the primary quarks and gluons. Furthermore, since the 
quarks and gluons carry non-zero color charges and the final hadrons do not, there can be 
no unique association of a jet of hadrons with a single initial quark or gluon. Nevertheless, 
with a suitable definition of the jet cross section one hopes to minimize the effect of long 
distance physics and of the inherent jet ambiguities and obtain a fairly precise picture of the 
short distance dynamics. 

Although the basic hard scattering processes studied in hadron-hadron and in electron- 
positron collisions are much the same, the overall event structure is quite different. In 
the e"^e~ case the initial state is purely electromagnetic and the entire final state can be 
thought of as arising from the short distance interaction of the virtual photon with the 
quarks. In this sense all of the hadrons in the final state are associated with the hard 
scattering process. In contrast the overall structure of the hadron-hadron case is much 
more complex. Of the large number of initial state partons, only one "active parton" from 
each incident hadron participates in the hard scattering process. Thus only a fraction 
of the hadrons in the final state are to be (loosely) associated with the hard scattering 
process, with the remainder ascribed to the "underlying event." This second contribution 
corresponds to the soft interactions of the remaining partons in the incident hadrons and, 
in first approximation, can be treated as uncorrelated with the hard process. The active 
partons also produce additional radiation in the form of initial state bremsstrahlung that is 
not present in e"'"e~ events. The underlying event plus the initial state radiation produce 
the characteristic "beam jets" of hadron collisions: particles with small momenta transverse 
to the beam axis, but possibly large momenta along the beam axis. The long distance 
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soft interactions responsible for the observed color singlet hadrons will, of course, result in 
some degree of dynamical coupling between all of these components. There will also be 
an essentially kinematical correlation induced by the fact that the jet selection or trigger 
process will generally be biased to choose events where the beam jets have higher than 
average global Et and multiplicity (i.e., the underlying event is noisier than average). 

These differences between the event structure of e'^e~ collisions and hadron-hadron col- 
lisions have, quite naturally, led to differences in the way jet definitions have been employed 
in the two cases. One might categorize the differences as follows. 

First, the cross sections studied are different. In e'^e~ collisions, where the entire event 
arises from an initially small number of energetic partons, one typically works with exclusive 
jet cross sections describing the production of exactly n jets and nothing else. In hadron- 
hadron collisions the practice has been to measure inclusive large pt jet cross sections, 
that is, cross sections to make n jets with specified properties plus any number of other 
unobserved jets or particles not in jets. 

Second, the variables used are different. For e+e" annihilation, one wants to emphasize 
rotational invariance. Thus the natural variables are energies E and polar angles 9, 0. 
For hadron-hadron collisions, one wants to emphasize invariance under boosts along the 
beam axis since, in fact, the cm. frame of the hard scattering is typically moving in the 
hadron-hadron cm. frame. Thus the natural variables are transverse momenta px or the 
corresponding "transverse energy" Et = £^sin^, azimuthal angle and pseudo-rapidity 
7] = -ln(tan(^/2)). 

Third, the jet definitions or algorithms employed to precisely define the (otherwise am- 
biguous) jets tend to be correspondingly different. As implied above in e'^e~ collisions one 
normally uses a jet definition that associates every final-state hadron uniquely with one of 
the jets. In hadron-hadron collisions producing high pt jets, only a small fraction of the 
final state hadrons are associated with the high pr jets. The other particles present in the 
event can be thought of as associated with the "beam jets" introduced earlier. One wants 
to keep the high pt jets distinct from the hadronic debris in the beam jets. For this reason, 
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one has typically used a cone definition [jT|, which was, in fact, inspired by the original the- 
oretical definition for jets in e^e~ collisions |0. A jet in this definition is a set of particles 
whose momentum vectors lie within a certain angular cone. Such a definition suppresses the 
effect of the beam jets, since only a small fraction of the low px particles in the beam jets 
will fall into the cone of a high pt jet. Furthermore, since this contribution is essentially 
determined by geometry, it is relatively easy to estimate. In the case of e+e^ collisions, the 
jet algorithm typically used experimentally is rather of the "successive combination" variety 
first introduced by the Jade group at DESY [^. In this kind of definition, one recursively 
groups sets of particles with "nearby" momenta, as defined by some measure, into larger 
sets of particles. The initial sets consist of just one particle each. The final sets are the jets. 

In this paper, we discuss a jet definition for hadron collisions that makes use of some 
of the ideas used in jet definitions in e^e'^ collisions. The definition is essentially that 
proposed by S. Catani, Yu. Dokshitzer, M. Seymour, and B. R. Webber [||, adapted for 
the measurement of inclusive, rather than exclusive, jet cross sections. We first define the 
algorithm. Then we comment on some of its features. Finally, we provide some evidence 
that this definition may have advantages compared to the cone definition that is currently 
the standard for hadron collisions. 

II. THE ALGORITHM 

We consider hadron collisions in the hadron-hadron cm. frame with the z-axis taken in 
the beam direction. We represent the final state of the collision as consisting of a starting 
set of "protojets" i with momenta . The starting p^ may be the momenta of individual 
particles, or each may be the total momentum of the particles whose paths are contained 
in a small cell of solid angle about the interaction point, as recorded in individual towers of 
a hadron calorimeter. In either case, we have in mind that the masses [PiPi^'^^'^ are small 
compared to the transverse momenta Pi^r-, so that the are essentially lightlike. Each 
protojet is characterized by its azimuthal angle 0j, its pseudorapidity rji = — ln(tan(6'j/2)), 
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and its transverse energy Et^i — \pT,i\- 

Starting with the initial hst of protojets, the jet algorithm recursively groups pairs of 
protojets together to form new protojets. The idea is that protojets with nearly parallel 
momenta should be joined, so that they will eventually form part of the same jet. The 
algorithm also determines when, for a particular protojet, joining should cease. This protojet 
is then labeled a completed "jet" and is not manipulated further. 

The algorithm depends on a parameter R, which should be chosen to be of order 1. This 
parameter is analogous to the cone size parameter in the cone algorithm. 

The algorithm begins with a list of protojets as described above and an empty list of 
completed jets. It then proceeds recursively as follows: 

1. For each protojet, define 



(1) 



and for each pair of protojets define 



(2) 



2. Find the smallest of all the di and dij and label it d, 



3. If dmin is a dij, merge protojets i and j into a new protojet k with: 



ET,k — ET,i + Etj 



(3) 



and 



Vk = [ET,iVi + EtjVj]/ ET,k, 



(4a) 



0fe = [ET,i(l>i + ET,j(t)j]/ET,k- 



(4b) 



4. If (imin is a di, the corresponding protojet i is "not mergable." Remove it from the list 
of protojets and add it to the hst of jets. 

5 



5. Go to step 1. 

The procedure continues until there are no more protojets. As it proceeds, it produces 
a hst of jets with successively larger values of di — E'^^^. 

III. COMMENTS 

The algorithm above produces a list containing many jets for each event. However, 
only the jets with large values of Et (which are the last to be added to the jet hst) are of 
much physical interest. The jets with smaller Et are "minijets" or just random debris from 
the beam jets. This situation is fine for the construction of an inclusive jet cross section. 
Consider, for instance, the one-jet inclusive cross section da/dEx for, say, Et — 100 GeV 
at ^/s = 1800 GeV. We first note that the high value of Et tells us that this jet is a signal 
of a short distance process. Second, we recall that in hadron collisions the probability for 
a parton collision decreases very quickly as the value of the parton-parton cm. energy \/§ 
increases. Thus it is very unlikely that an event with a 100 GeV jet contains other high 
Et jets beyond a second jet with Et ~ 100 GeV that is needed to balance the transverse 
momentum. That is to say, the one jet inclusive cross section is primarily sensitive to the 
highest Et jets in hard scattering events, even though the jet list for each event contains 
many jets covering a wide range of transverse energies. 

It is crucial that a jet cross section be "infrared safe." That is, when the cross section is 
calculated in QCD at the parton level, the cross section must be finite, despite the infrared 
divergences present in the Feynman diagrams. At the level of the physical hadrons, this 
means that the cross section is not sensitive to long distance effects. The infrared divergences 
in Feynman diagrams come from configurations in which a parton emits a soft gluon, with 
q'^ 0, or in which an outgoing parton divides into two coUinear partons, or in which 
an incoming parton emits another parton that carries away a fraction of its longitudinal 
momentum but no transverse momentum. The probability for one of these configurations to 
occur is infrared sensitive, and infinite in fixed order perturbation theory. However, unitarity 
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dictates that the sum of the probabihties for one of these configurations to happen or not to 
happen is 1. For this reason, infrared safety is achieved if the measured jet variables do not 
change when an Et —>■ parton is emitted or when a parton divides into colhnear partons. 
We note that da/dEx and similar jet cross sections produced using the algorithm described 
in the previous section will have this property. If a parton divides into two partons with 
collinear momenta, then the algorithm immediately recombines them, producing the same 
result as if the parton had not divided. Similarly, an Et parton may wind up in one of 
the high Et jets or it may be left by itself, but in the limit that its Et tends to zero it does 
not affect the transverse energy or angle of the high Et jet. 

We have discussed the one jet inclusive cross section. For the two jet inclusive cross 
section, it would be sensible to pick the two jets in each event that have the highest transverse 
energies, just as has been done for jets defined with the cone algorithm. By analogy, in 
defining the one jet inclusive cross section one might be tempted to pick only the most 
energetic jet in each event. However, at the Born level there are two jets with exactly equal 
transverse energies. Which one winds up with the most transverse energy is affected by a 
long-distance process, the emission of very low Et gluons. Thus the resulting cross section 
would not be infrared safe. 

The algorithm defined above is very similar to the simplest version of the various options 
discussed by Catani et al. 0], which are generalizations of the "Durham" algorithm for 
e^e~ annihilation. The differences arise primarily from questions of emphasis. Catani et 
al. take the approach that the analysis should be kept as similar as possible to earlier e^e~ 
work. Thus they fix the parameter R at the value 1 and focus on exclusive ]ei cross sections. 
They also introduce two further parameters. The first additional parameter is dent- When 
the smallest di or dij is larger than dent-, their recursion halts. The jets that have been 
generated thus far, all of which have E'^^ < dcuti are regarded as part of the beam jet. The 
remaining protojets, all of which have Eti > (icut, are treated as resulting from the hard 
scattering process. These protojets are then resolved into the final "jets" in direct analogy to 
the e+e~ case using a resolution parameter ?/cut- This is a sensible way to define an exclusive 
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jet cross section in analogy to the e~^e~ case but now including the beam jets. For instance, 
one might measure in this way a cross section cr3((icut; Z/cut) to produce exactly three high 
Et jets plus the beam jets. However, we prefer to maintain a similarity with the previous 
cone algorithm work in hadron-hadron collisions. The jet definition in the preceding section 
is intended to define inclusive jet cross sections in terms of the single angular resolution 
parameter R (which plays a role similar to ?/cut)- For example, the two jet inclusive cross 
section da/dMjj defined in this way is a function of only the jet-jet invariant mass Mjj and 
R. It would be an additional complication if it also depended on the parameters dcut and 

I/cut- 

We note, however, that an exclusive n-jet cross section can be defined using the algorithm 
in Section 2. In this case, one needs a jet hardness parameter to play the role of dent- A 
convenient choice is to count only jets with transverse energies above a cutoff Extent- Except 
for the issue of variable R, this is essentially the i/cnt = 1 scenario of Catani et al. 

The function dij given in eq. measures how "nearby" the pair of protojets {i,j) is. 
As in other algorithms of the successive combination type, the idea is to combine first the 
protojets that are "nearest", and thus have the smallest dij. There are, of course, other 
possibilities for the function dij, the measure of "nearness". For instance, one might use the 
invariant pair mass [for massless protojets i and j], 

<5 = M^j = 2ET,iET,i[cosh{r], - r],) - cos(0i - 0,)]. (5) 

For small \rii — rjjl and |0j — (f)j\ this is 

^ ET,ErAiv^ - + (0. - <P,n (6) 

Such a choice for dij yields an algorithm analogous to the "Jade" version of the successive 
combination algorithm used in e+e^ annihilation. The corresponding algorithm with the 
factor Ex^iExj replaced by m.m{E^^^, E^ j) , as in eq. (0), is analogous to the "Durham" 
algorithm |^ in e~^e~ annihilation as noted earlier. A discussion of the relative merits of 
these choices in the context of e^e~ annihilation may be found in ||^. The Durham algorithm 
has been also discussed in the context of lepton-hadron collisions in . 
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IV. COMPARISON WITH CONE ALGORITHM 



Now consider the algorithm advocated here. At any stage in the operation of the algo- 
rithm, the two protojets i and j with the smallest value of dij are merged if dij of eq. (0) is 
less than the smaller of E^^ and E^j. That is, they are merged if 

^/(r/, - Vj)' + (0. - < R- (7) 

Thus the issue of which protojets to merge first depends on transverse energies and angles, 
but the issue of whether to merge two protojets or to declare that they cannot be merged is 
solely a question of the angle between them. 

The merging condition in eq. makes it clear that the present algorithm is in fact not 
so different from a cone algorithm. The latter is typically defined in terms of the particles 
n whose momenta lie within a cone centered on the jet axis {rjj, (pj) in pseudorapidity rj 
and azimuthal angle 0, 

^{Vn-Vjy + i<Pn-M' <R- (8) 

The jet angles {rjj, are the averages of the particles' angles, 

Vj= J2 PT,nVn/ET,j, (9a) 

nGcone 

0J= PT,n(j)n/ET,J , (9b) 

nGcone 

with 

Et,J= E PT,n, (10) 
n Scone 

analogous to eqs. (^ and (0), respectively. This process is iterated until the cone center 
matches the jet center (?7j,0j) computed in eq. (^. 

The definition of the jet axis given in eqs. (^ and (|T0|) is chosen to be simple when 
expressed in the natural variables {E^, r], 0). Of course, other choices could also yield infrared 
safe jet definitions. Thus this definition should be regarded as a convenient convention. This 
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convention has been continued in the merging conditions, eqs. (|^) and (|), of the successive 
combination algorithm. 

While the successive combination algorithm never assigns a particle to more than one 
jet, this is not the case for the cone algorithm as defined so far. It is possible for jet cones 
to overlap, so that one particle is contained in more than one jet. This issue was discussed 
in P] in the context of the order perturbative calculation. At this order it is possible 
for two one-parton jets to lie within the cone of a two parton jet. In the calculation [§, 
such jets are "merged." That is, the two parton jet is kept and the one parton subjets are 
not considered as being legitimate jets on their own. In a physical event, with many more 
particles, the merging question is more serious, and a criterion for merging must be part of 
the experimental algorithm [p!0|j . 

The difference between cone and successive combination algorithm jets is apparent even 
in the simplest example of merging two partons (or hadrons) to make a jet. In the cone 
algorithm [Q, two partons i and j are merged if each falls within an angular distance R of 
the jet axis defined by eq. (^^. The parton with the smaller E^, call it is farther from 
the jet axis. In the cone algorithm, the limit on angular separation between this parton and 
the jet J (including the partons i and j) has the simple form 



The corresponding relation in terms of the angular separation of the two partons is 



The limit on the right hand side lies in the range R < {Etj/[Et,j — ET,i]}R < 2R since 
< Ex^i < Et^j/2. Thus configurations are possible with two equally energetic partons 
located near opposite edges of the cone, and nothing in the center of the cone. These are 
precisely the configurations where the merging issue arises in order perturbation theory. 

In the successive combination algorithm, it is the separation between the two partons 
(or protojets) that has the simple limit of eq. (0). It is this angle between i and j, not 




(11) 




(12) 
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that between i and the final jet direction, that controls the question of merging. The 
corresponding relation for the angular separation between parton i (the lower E-r parton), 
and the jet J in order for merging to occur is 



^J{v^-VJy + {<P^-<pJy< ^'''' ^^' R, (13) 



where the right hand side is in the range R/2 < [Et,j — ET,i\/ET,j < R- Thus the lower-i?T 
parton can be far from the jet axis, up to a maximum separation R, while the higher Et 
parton must be closer to the jet axis. 

Because of the difference between eq. (|ll]) and eq. (|13D, the average distribution of 
transverse energy within jets depends on which algorithm one uses. With the successive 
combination algorithm, there is less transverse energy near the edge of the allowed angular 
region than there is with the cone algorithm. This is illustrated in Fig. 1. The quantity 
plotted is the order perturbative result for the average transverse energy fraction as a 
function of distance from the center of the jet, for jets with Et = 100 GeV at ^/s = 1800 
GeV for i? = 1. The histogram represents the Et fraction in angular annuli, r to r + 0.1, 



where r is the distance from the jet center, r = y (r/ — rjjy + (0 — 0j)^. Thus the sum over 
all bins, r < R including < r < 0.1, should equal unity. Note also that the calculation 
yields energy outside of the cone, r > R. This arises in the perturbative calculation from 
configurations with only one parton in the cone (the jet) but with another parton nearby. 

The cone algorithm, in principle, requires one to find all possible solutions of the cone 
matching conditions, eqs. (H), (^ and (p!0D , before beginning the merging algorithm. In 
experimental applications [|1^, however, one may begin the search for valid jet cones using 



"initiator" calorimeter cells, i.e., cells with Et above some threshold value. In this case, 
possible jets that consist of two widely separated subjets may not be recognized because 
there is no initiator cell between the subjets. Thus, the two widely separated jets may not 
be merged. The theoretical study suggests that the jet finding algorithm used by the 
CDF group is likely not to merge two subjets when their separation is greater than a value 



Rsep, which appears to be about 1.3R in practice as is consistent with CDF studies [|10 
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In this language, the successive combination algorithm, when applied to two parton jets, 
corresponds to the hmiting case Rgep = R- 

At higher orders in perturbation theory or for more realistic jets containing many par- 
ticles, it is eq. ( pA]) that survives in the cone algorithm case and eq. (|^) for the successive 
combination algorithm. Thus cone jets always have well defined, smooth boundaries, al- 
though the amount of Et near the edge will depend in practice on how the issue of the 
merging of overlapping jets is handled. For the successive combination algorithm there is 
no "merging question." It is automatically dealt with by the algorithm as all particles are 
assigned to a unique jet. Furthermore, only small Et particles can be as far as R from the 
jet axis. However, there is a price for this simplicity. Particles with very small E^ can be 
very far from the jet axis (> R). Thus jet boundaries can be complicated. Here, the issue 
is intertwined jets rather than overlapping jets. 

Does it follow that the successive combination algorithm is better in some sense than the 
cone algorithm, or is it just different? The criteria that we will examine are 1) estimated size 
of higher order perturbative corrections to a jet cross section calculated using fixed order 
perturbation theory, 2) estimated size of corrections to the calculated cross section arising 
from "power suppressed" or "hadronization" effects, and 3) simplicity and definiteness of 
the algorithm. 

We turn first to the estimated size of perturbative corrections. We consider as an example 
the one jet inclusive cross section da/dETdr] at -^i = 1800 GeV, averaged over the rapidity 
range 0.1 < |?7| < 0.7 used by the CDF Collaboration where the comparison of theory 



with data has been very successful [§,|l3]. Again we calculate this cross section at order 



according to both jet definitions, using the methods described in ||^ and 0. The question is 
then, how big are the order corrections? Of course, we can only make estimates. One way 
to do this is to examine the dependence of the calculated cross section on the renormalization 
and factorization scale fi. If we could calculate to order a^, then the renormalization group 
guarantees that some of the terms would cancel most of the /i dependence that occurs 
in the cross section (so that the derivative of the cross section with respect to ^ would 
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then be of order a^). Thus the contributions to the cross section are hkely to be at least 
as large as the difference between the cross section calculated with a "best guess" for /i, 
which we take to be = Et/2, and the cross section calculated with the alternative choices 
fi = Et/4: or yU = Et- In Fig. 2a, we show the cross section at Et = 100 GeV as a function 
of the clustering parameter R calculated according to both algorithms for these three values 
of /i. 

This rather busy graph becomes remarkably simple, as in Fig. 2b, if we rescale the 
parameter R in the two cases. Let the clustering parameter in the successive combination 
algorithm, eq. (H), be called -Rcomb for the present purposes and plot the cross section versus 

R' = -Rcomb- (14) 

For comparison, we label the clustering parameter in the case of the cone algorithm, eq. (^, 
as -Rcone and in Fig. 2b we display the cone results as a function of a scaled parameter 

R' = 1.35 X i?,o,e. (15) 

We see from the graph that, for each value of /x, the calculated cross sections are almost 
identical as long as we identify -Rcomb ~ 1-35 x -Rcone- Thus, for instance, a jet cross section 
calculated or measured with the cone algorithm using the standard value -Rcone = 0.7 should 
be compared to a jet cross section with the successive combination algorithm using -Rcomb = 
0.945 ^ 1.0. 

In Fig. 3 we plot the order inclusive jet cross section as defined by the successive 
combination algorithm with -Rcomb = 1-0 versus the jet E^ at y/s = 1800 GeV taking /i = 
Et/2. The corresponding cross section calculated using the cone algorithm with -Rcone = 0.7 



is nearly identical (see, for example. Fig. 1 of |T2[)) and is not shown as it would not be 
distinguishable. In the range Et/4: < jj, < E^, 10 GeV < E^ < 500 GeV, we find that the 
two algorithms give order al theoretical cross sections that agree to within 10%. 

We conclude from this comparison that the successive combination algorithm is neither 
better or worse than the cone algorithm, at least as judged according to this /i-dependence 
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standard. We also note that the argument that -Rcone ~ 0.7 is a particularly stable and 
thus sensible regime for comparing fixed order perturbation theory with experiment is now 
translated into -Rcomb ~ 1-0. This is just the regime studied by Catani et al. [Q. 

We now examine the question of higher order corrections from a point of view that 
relates directly to the reasoning behind the successive combination algorithms: the idea of 
putting parton showers back together to reconstruct the parent parton. We suggest that jet 
cross sections defined with the cone algorithm may have larger higher order perturbative 
corrections than do jet cross sections defined with the successive combination algorithm 
because of what may be called "edge of the cone" effects. Consider the application of 
the cone algorithm to two partons, 1 and 2, with roughly equal transverse energies and 
an angular separation of approximately 2_Rcone, the troublesome configuration mentioned 
earlier. These two partons are near the edge of the region in which they can form an 
allowable cone jet. Suppose that parton 2 splits into two partons, 2a and 2b, that each have 
substantial transverse energy. If the angle separating partons 2a and 2b is infinitesimal, 
then they will both fit into the jet cone. The resulting jet (l,2a,2b) will have the same 
direction and transverse energy as the jet (1,2) that one obtains if the splitting did not 
occur. However, if the angle separating partons 2a and 2b is small but not infinitesimal, 
it can very well be that the three partons (1,2a, 2b) cannot fit into a cone to form a single 
jet. Since the matrix element for such a parton splitting is large (although not infinite), one 
may worry that there will be corresponding large order corrections to jet cross sections 
calculated with the cone algorithm. What is the situation with the successive combination 
algorithm? Here we should consider partons 1 and 2 separated by an angle of approximately 
-Rcomb, near the edge of the region in which they can form an allowable jet. If parton 2 splits 
into partons 2a and 2b with a small angular separation, then the successive combination 
algorithm will combine them together into a protojet 2' that approximates parton 2. As 
long as the angle separating partons 2a and 2b is not too large, the protojet 2' will have jet 
parameters close enough to those of parton 2 that the algorithm will then combine 2' with 
parton 1. On the basis of this argument, we expect that order corrections to jet cross 
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sections calculated with the successive combination algorithm may be smaller than with the 
cone algorithm. For similar reasons, we expect also that jet cross sections calculated with 
the successive combination algorithm will also exhibit smaller corrections attributable to the 
final combination of partons into hadrons. Unfortunately, from the order perturbative 
calculation we cannot determine the magnitude of these edge of the cone effects. [| 

Finally, we comment on the relative merits of the two algorithms from the point of view of 
simplicity and definiteness. Here the cone algorithm appears at first to have the advantage. 
With the cone algorithm, a jet consists simply of all the particles whose momentum vectors 
fit into an apparently well defined and regular cone centered on the jet axis. This is a simple 
and appealing idea. The successive combination algorithm, in contrast, takes more effort to 
define and does not yield regular jet shapes in the f] — (p plane. However, in the application 
of the cone algorithm, one quickly discovers that ambiguous cases with overlapping jet cones 
arise. What do you do when two or more jets contain particles in common? The algorithm 
must be expanded to cover these cases. Unfortunately, there are many ways to proceed, and 
none of them is particularly simple. Thus the cone algorithm actually consists of a simple 
part that works beautifully for joining two partons into a jet and a complicated part that 
one is forced to use in order to deal with real-world multi-hadron events. With the successive 
combination algorithm, the recursive steps require some thought to define, but once they 
are set, the whole algorithm is complete. 

V. SUMMARY 

We have discussed a jet definition for inclusive jet measurements in hadron collisions that 
makes use of ideas previously applied to jet definitions for e"*"e~ collisions. The algorithm 

^This is clearly a subject for Monte Carlo study as in Q. Unfortunately the cone algorithm used 
in that analysis is quite different from that described here. In particular, the cone algorithm in |Q] 
is not infrared safe so that the results presented are difficult to interpret. 
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identifies a jet by successively combining "nearby" pairs of particles or protojets. The 
concept of "nearby" is measured in {Et, rj, 0) space and involves a limit i?comb on angular 
separations that is very similar to the usual cone algorithm parameter Rcone- We find 
that, calculated perturbatively, the inclusive jet cross section that results from the new 
algorithm with parameter i?comb is essentially identical to the cone algorithm result with 
-Rcone = -Rcomb/1-35. While the final geometry of a jet defined by the successive combination 
algorithm is likely to be more complex than that from the cone algorithm, the former 
definition has the advantage that there is no problem with overlapping jets as there is 
in the cone case. We have also presented a qualitative argument that, due to "edge of 
the cone" effects, cross sections calculated with the successive combination algorithm are 
likely to exhibit smaller higher order and hadronization corrections. Only further detailed 
experimental and theoretical studies can demonstrate whether the successive combination 
type algorithm has quantitative advantages over other algorithms. 
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FIGURES 

FIG. 1. Fraction of jet Et in angular annuli r to r + 0.1 comparing the cone algorithm with the 
the successive combination case. In both cases the jet has R = 1.0, Ex = 100 GeV, ^/s = 1800 GeV, 
0.1 < \r]j\ < 0.7 with renormalization/factorization scale fi = Et/2 and the structure functions of 
HMRS(B) |ll|. 

FIG. 2. Order inclusive jet cross section for Et = 100 GeV, ^/s = 1800 GeV, averaged 
over rjj in the range 0.1 < |r/j| < 0.7 with the structure functions of HMRS(B) for the two 
algorithms specified in the text. The curves for the cone algorithm are: /j, = Et (solid), /i = Et/2 
(dot-dash), = Et/^ (dot-dot-dot-dash); for the successive combination algorithm: ^ = Et 
(long dash), /x = Et/2 (medium dash), /x = Et/^ (short dash), a) Standard case plotted versus 
R = -Rconc = -Rcomb- b) Same as a) except that the cone algorithm is plotted versus R' = 1.35i?cone 
while the successive combination case has R' = Rcomb- 

FIG. 3. Order inclusive jet cross section as defined by the successive combination algorithm 
with -Rcomb = 1-0 versus the jet Et for ^/s = 1800 GeV, = Et/^, averaged over rjj in the range 
0.1 < \r]j\ < 0.7 with the structure functions of HMRS(B) |ll|. 
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